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1. Introduction 

It is not yet feasible to perform simulations on lattices that can simultaneously represent the 
two relevant scales of B-physics: the low energy scale Aqcd> requiring large physical lattice size, 
and the high energy scale of the b-quark mass nib, requiring very small lattice spacing a. 

A promising alternative is to consider (lattice) heavy-quark effective theory (HQET), which 
allows for an elegant theoretical treatment, with the possibility of fully nonperturbative renormal- 
ization [1] (see [2] for a review). The approach is briefly described as follows. HQET provides 
a valid low-momentum description for systems with one heavy quark, with manifest heavy-quark 
symmetry in the limit nib — ► °°- The heavy-quark flavor and spin symmetries are broken at fi- 
nite values of rtib respectively by kinetic and spin terms, with first-order corrections to the static 
Lagrangian parametrized by (O^m and G) sp i n 

JZ? H Q ET = Jjf h ( x ) D y h (x)- (Okin ^ kin - W spin ^ spin , (1.1) 



where 



? kin = Y h (x)D 2 \jf h (x) , ^ spin = \j/ h (x)a-BYh{x). (1.2) 




These 0(\/nib) corrections are incorporated by an expansion of the statistical weight in l/m& 
such that i^ki n , ^spin are treated as insertions into static correlation functions. This guarantees the 
existence of a continuum limit, with results that are independent of the regularization, provided that 
the renormalization be done nonperturbatively. 

As a consequence, expansions for masses and decay constants are given respectively by 

m B = m bare + £ stat + a^E** + w spin £ spin (1.3) 

and 

= zf2 E V at (l + c* QE V A + «W kin + «s P inP spin ) » d-4) 

where the parameters nibare and Z^ ET are written as sums of a static and an 0(1 /nib) term (denoted 
respectively with the superscripts "stat" and "1 /nib" below), and c^ QET is of order l/m# (see e.g. [3] 
for unexplained notation). Bare energies (E stat , etc.) and matrix elements (p stat , etc.) are computed 
in the numerical simulation. 

The divergences (with inverse powers of a) in the above parameters are cancelled through the 
nonperturbative renormalization, which is based on a matching of HQET parameters to QCD on 
lattices of small physical volume — where fine lattice spacings can be considered — and extrapo- 
lation to a large volume by the step-scaling method. Such an analysis has been recently completed 
for the quenched case [4]. In particular, there are nonperturbative (quenched) determinations of 
the static coefficients n^are anc ^ ^A tat f° r HYP1 and HYP2 static-quark actions [5] at the physical 
b-quark mass, and similarly for the 0(1 /nib) parameters (O^n, ft) sp j n , ^bare*' z\' mb and c^ ET . 

The newly determined HQET parameters are very precise (with errors of a couple of a percent 
in the static case) and show the expected behavior with a. They are used in our calculations reported 
here, to perform the nonperturbative renormalization of the (bare) observables computed in the 
simulation. Of course, in order to keep a high precision, also these bare quantities have to be 
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accurately determined. This is accomplished by an efficient use of the generalized eigenvalue 
problem (GEVP) for extracting energy levels E„ and matrix elements, as described below. 

A significant source of systematic errors in the determination of energy levels in lattice simu- 
lations is the contamination from excited states in the time correlators 

oo 

c(t) = (o(t)o(o)) = £ |Hd|o>| 2 ^ (i.5) 

n=\ 

of fields 0(t) with the quantum numbers of a given bound state. 

Instead of starting from simple local fields O and getting the (ground- state) energy from an 
effective-mass plateau in C(t) as defined above, it is then advantageous to consider all-to-all prop- 
agators [6] and to solve, instead, the GEVP 

C(t)v n (t,t ) = K(t,to)C(to)v n (t,t ), (1.6) 

where t > to and C(t) is now a matrix of correlators, given by 

oo 

Qj{t) = (Oi(t)Oj(0)) = ^e-^M^-, i,j = \,...,N. (1.7) 

n=l 

The chosen interpolators Oj are taken (hopefully) linearly independent, e.g. they may be built from 
the smeared quark fields using N different smearing levels. The matrix elements VP n / are denned by 

^ni = {^n)i = (n\6i\0) , (m\n) = 8 mn . (1.8) 

One thus computes Qj for the interpolator basis O, from the numerical simulation, then gets 
effective energy levels Ef^ and estimates for the matrix elements *P„,- from the solution X n (t,to) of 
the GEVP at large t. For the energies 

£fM) = -lo gj M^ (L9) 
a A„(t + a,to) 

it is shown [7] that £^ ff (?,?o) converges exponentially as t —> oo (and fixed to) to the true energy E n . 
However, since the exponential falloff of higher contributions may be slow, it is also essential to 
study the convergence as a function of to in order to achieve the required efficiency for the method. 
This has been done in [3] , by explicit application of (ordinary) perturbation theory to a hypothetical 
truncated problem where only ,/V levels contribute. The solution in this case is exactly given by the 
true energies, and corrections due to the higher states are treated perturbatively. We get 

Ef(t,t ) = E n + e n (t, to) (1.10) 

for the energies and 

oo 

e-*(Jf (Mo)) + |0) = |«>+£ *w(f,fo) l«> d- 11 ) 

n'=\ 

for the eigenstates of the Hamiltonian, which may be estimated through 1 

Jf(Mo) = R n (0,v n (t,to)), (1.12) 

X n (to + a,to) 



. i £/2 

*n = (Vn(tM,C(t)Vn(t,to))- / 



(1.13) 



'The choice of R„ we make here leads to smaller statistical errors than the one in [3], while the form of the correc- 
tions K remains unchanged. 
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In our analysis we see that, due to cancellations of ?-independent terms in the effective energy, 
the first-order corrections in e n (t,to) are independent of to and very strongly suppressed at large 
t. We identify two regimes: 1) for to < t/2, the 2nd-order corrections dominate and their expo- 
nential suppression is given by the smallest energy gap \E m — E n \ = AE„ hn between level n and its 
neighboring levels m; and 2) for to > t/2, the lst-order corrections dominate and the suppression is 
given by the large gap AEn+i , n - Amplitudes K nn i(t,to) get main contributions from the first-order 
corrections. For fixed t — to these are also suppressed with AZsjv+i , n - Clearly, the appearance of 
large energy gaps in the second regime improves convergence significantly. We therefore work 
with t, to combinations in this regime. 

A very important step of our approach is to realize that the same perturbative analysis may be 
applied to get the 1 /nib corrections in the HQET correlation functions mentioned previously 

Qj(t) = qf(t) + coC)' m \t) + 0(co 2 ) , (1.14) 

where the combined 0(1/ m&) corrections are symbolized by the expansion parameter CO. Following 
the same procedure as above, we get similar exponential suppressions (with the static energy gaps) 
for static and 0(1 /nib) terms in the effective theory. We arrive at 

Ef(t,t Q ) = Ef> sta \t,to) + (oEf' l/m "(t,t ) + 0((0 2 ) (1.15) 

with 

Ef' sm (t,t ) = E%* + |fe^' + . . . , (1.16) 

Ef' 1/m »(t,t ) = E l J mb + - ft*t A45>-^V + . . . . (1.17) 

and similarly for matrix elements. 

Preliminary results of our application of the methods described in this section are presented 
next. A more detailed version of this study will be presented elsewhere [8]. 

2. Results 

We carried out a study of static-light B s -mesons in quenched HQET with the nonperturbative 
parameters described in the previous section, employing the HYP1 and HYP2 lattice actions for 
the static quark and an 0(a) -improved Wilson action for the strange quark in the simulations. The 
lattices considered were of the form L 3 x 2L with periodic boundary conditions. We took L rs 1.5 
fm and lattice spacings 0.1 fm, 0.07 fm and 0.05 fm, corresponding respectively to j8 = 6.0219, 
6.2885 and 6.4956. We used all-to-all strange-quark propagators constructed from approximate low 
modes, with 100 configurations. Gauge links in interpolating fields were smeared with 3 iterations 
of (spatial) APE smearing, whereas Gaussian smearing (8 levels) was used for the strange-quark 
field. A simple 7075 structure in Dirac space was taken for all 8 interpolating fields. Also, the local 
field (no smearing) was included in order to compute the decay constant. 

The resulting (8 x 8) correlation matrix may be conveniently truncated to mN xN one and 
the GEVP solved for each N, so that results can be studied as a function of N. We have considered 
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Figure 1: Static energy spectrum obtained from the first basis of interpolators for the case of HYP2 action 
and j3 w 6.3: ground state (left) and first excited state (right). 
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Figure 2: Static energy spectrum obtained from the second basis of interpolators for the case of HYP2 
action and j3 « 6.3: ground state (left) and first excited state (right). 



two bases of interpolators. The first basis was obtained by projecting (at f; 
eigenvectors of C(?,) with the largest eigenvalues 



0.2 fm) with the N 



C(ti)b n = X n b n 



C (NxN) (t) 



b\C{t)b m , n,m<N. 



(2.1) 



For N not too large, this avoids numerical instabilities and large statistical errors in the GEVP. This 
basis was used in [3], where also the normalization of the different smeared fields is specified. 
Note that the (relative) normalization does matter in Eq. (2.1). The second basis was picked from 
unprotected interpolators, sampling the different smearing levels (from 1 to 7) as {1,7}, {1,4, 7}, 
etc. Our results for the effective energies using the two bases are shown in Figs. 1 and 2, where 
the solid lines correspond to a simultaneous fit of the various energy levels and values of N to the 
behavior 2 in Eq. (1.16). In both cases, we observe the predicted to independence in the GEVP 
solution for the energies. We see a much stronger dependence on N for the ground state in the first 
case. Nevertheless, our final results remain unchanged within their errors, which are determined as 
described in [3]. 

We then take the continuum limit, as shown in Figs. 3 and 4. We see that the correction to 
the ground-state energy due to terms of order 1 /mi,, which is positive for finite a, is quite small 



2 For N = 2 care must be taken, since the subleading corrections may not be negligible in the available / interval. 
This happens for our first basis above. In any case, we do not include data from N = 2 in our analysis. 
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Figure 3: Continuum extrapolation (joint limit for HYP1, HYP2) of the energy splittings (in MeV). Shown 
are AE21 (lower two curves, respectively the static and the full values) and AE31 (static, uppermost curve). 
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Figure 4: Continuum extrapolation (joint limit for HYP1, HYP2) of the pseudoscalar meson decay constant 
in the static limit (left) and to order 1 jm\, (right). 

(consistent with zero) in the continuum limit. Our results for the pseudoscalar meson decay con- 
stant, both in the static limit and including 0(1 /nib) corrections, are shown in terms of the com- 
bination <J> H( 2 ET = Fps ^/mps/Cps, where Cps(M/Aqcd) is a known matching function and <J> RGI 
denotes the renormalization-group-invariant matrix element of the static axial current [9]. These 
two continuum extrapolations are shown in comparison with fully relativistic heavy-light (around 
charm-strange) data from [9] in Fig. 5 below. Note that, up to perturbative corrections of order a 
in Cps, HQET predicts a behavior const. + 0(1 /romps) in this graph. Surprisingly no l/(romps) 2 
terms are visible, even with our rather small errors. 



3. Conclusions 



The combined use of nonperturbatively determined HQET parameters (in action and currents) 
and efficient GEVP allows us to reach precisions of a few percent in matrix elements and of a 
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Figure 5: Comparison of the continuum values for the pseudoscalar meson decay constant from Fig. 4 to 
fully relativistic data in the charm region. The solid line is a linear interpolation between the static limit and 
the points around the charm-quark mass, which corresponds to 1/romps « 0.2. 

few MeV in energy levels, even with only a moderate number of configurations. The method is 
robust with respect to the choice of interpolator basis. All parameters have been determined non- 
perturbatively and in particular power divergences are completely subtracted. We see that HQET 
plus 0(l/mfo) corrections at the b-quark mass agrees well with an interpolation between the static 
point and the charm region, indicating that linearity in l/m extends even to the charm point. A 
corresponding study for Nf = 2 is in progress. 
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